Astronomy & Astrophysics manuscript no. 
February 28, 2013 



18658 



©ESO2013 



(N 

o 

(N 



O 
(N 

O 



> 

\o 
in 
^, 

rn 

O 



EVN observations of 6.7 GHz methanol maser polarization in 

massive star-forming regions. 

G. Surcis *'''^, W.H.T. Vlemmings^'^, H.J. van Langevelde^'^, and B. Hutawarakom Kramer^'^ 

' Argelander-Institut fiir Astronomie der Universitat Bonn, Auf dem Hiigel 71, 53121 Bonn, Germany 

e-mail: gsurcis@astro.uni-bonn.de 

- Joint Institute for VLB! in Europe, Postbus 2, 7990 AA Dwingeloo, The Netherlands 

^ Chalmers University of Technology, Onsala Space Observatory, SE-439 92 Onsala, Sweden 

* Sterrewacht Leiden, Leiden University, Postbus 9513, 2300 RA Leiden, The Netherlands 
' Max-Planck Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 

* National Astronomical Research Institute of Thailand, Ministry of Science and Technology, Rama VI Rd., Bangkok 10400, Thailand 



Received ; accepted 



ABSTRACT 



Context. The role of magnetic fields in the formation of high-mass stars is still under debate, and recent measurements of their 
orientation and strength by using polarized maser emissions are contributing new insights. Masers polarization, in particular of the 
6.7-GHz methanol masers, are one of the best probes of the magnetic field morphologies around massive protostars. 
Aims. Determining the magnetic field morphology around an increasing number of massive protostars at milliarcsecond resolution 
by observing 6.7-GHz methanol masers is crucial to better understand the role of magnetic fields in massive star formation. 
Methods. The First EVN Group consists of 4 massive star-forming complexes: W51, W48, IRAS 18556H-0138, and W3(0H). These 
contain well-studied Hii regions from some of which molecular bipolar outflows were also detected (W51-e2, G35.20-0.74N). Nine 
of the European VLBI Network antennas were used to measure the linear polarization and Zeeman-splitting of the 6.7-GHz methanol 
masers in the star-forming regions of the First EVN Group. 

Results. We detected a total of 154 CH3OH masers, one third of these towards W3(0H). Fractional linear polarization (1.2-1 1.5%) 
was detected towards 55 CH3OH masers. The linear polarization vectors are well-ordered in all the massive star-foiTning regions. We 
measured significant Zeeman-splitting in 3 massive star-forming regions (W51, W48, and W3(OH)) revealing a range of separations 
-3.5 ms"'< AVz < 3.8 ms"' with the smallest lAV^I = 0.4 ms"'. We were also able to compare our magnetic field results with those 
obtained from submillimeter wavelength dust observation in W51 and show that the magnetic field at low and high resolutions are in 
perfect agreement. 
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1. Introduction 

. 5^ A number of different scenarios have been proposed to explain 
k^ the formation of stars with masses larger than about 8 Mq. 
;_] , These include formation through the merger of less massive 
Ci • stars (Coalescence model; e.g.. Bally & Zinnecker 2005) or 
through the accretion of unbound gas from the molecular cloud 
(Competitive accretion model; e.g., Bonnell et al. 2004). In the 
third scenario, the Core accretion model, massive stars form 
through gravitational collapse, which involves disc-assisted ac- 
cretion to overcome radiation pressure (e.g., McKee & Tan 
2003). This scenario is similar to the favored picture of low -mass 
star formation, in which magnetic fields are thought to regulate 
the collapse, to transfer the angular momentum and to power 
the bipolar outflows. Moreover, during the formation of low- 
mass stars the gravitational collapse of molecular clouds pro- 
ceeds preferentially along the magnetic field lines, giving rise 
to large rotating disc or torus structures orthogonal to the mag- 
netic field (e.g., Matsumoto & Tomisaka 2004). Consequently, 
the molecular bipolar outflows, which originate from the proto- 
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star, are driven parallel to the original magnetic field orientation. 
Although the magnetic fields play such an important and cru- 
cial role in the formation of low-mass stars, their role in the for- 
mation of high-mass stars is still under debate. Fairly recently 
magnetic fields have been included in the theoretical simula- 
tions of high-mass star formation (e.g. Banerjee & Pudritz 2007). 
Therefore, providing new measurements of magnetic fields ori- 
entation and strength at milliarcsecond (mas) resolution is fun- 
damental to understand the formation process of high-mass stars. 

Current dust observations of magnetic fields in massive star- 
forming regions are often limited to low density regions and/or 
envelopes at scales of several thousands astronomical units (e.g., 
Koch et al. 2010). Linear polarization observations of dust also 
only provide information on the magnetic field in the plane of 
the sky, so these observations have been yet unable to probe the 
strength and the full structure of the magnetic field close to the 
protostars and around protostellar discs (e.g.. Tang et al. 2010). 

The best probes of magnetic fields in the high density regions 
close to massive protostars currently available are masers . Their 
bright and narrow spectral line emission is ideal for measuring 
the Zeeman-splitting even though the exact proportionality be- 
tween the measured splitting and the magnetic field strength is 
still uncertain (see Vlemmings et al. 201 1), as well as for deter- 
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mining the orientation of the magnetic field in 3-dimension (e.g., 
Vlemmings et al. 2010). Recently 6.7-GHz CH3OH masers, 
which are the most abundant maser species in massive star- 
forming regions, have started to play a crucial role in determin- 
ing the magnetic field morphology close to massive protostars 
(e.g., Vlemmings et al. 2006, 2010, Surcis et al. 2009, 2011a). 
The CH3OH masers probe magnetic fields along outflows and on 
surfaces of tori (Surcis et al. 2009;Vlemmings et al. 2010; Surcis 
et al. 2011a). Therefore, determining the magnetic field mor- 
phology around a large number massive protostars by observing 
the polarized emission of 6.7-GHz CH3OH masers is of funda- 
mental importance. Here we show the results of our First EVN 
Group composed of 4 massive star-forming regions, which are 
described in details in Sect. 2. Observations of a second group 
of 5 sources are in the process of being analyzed. 



2. The First EVN Group 

We selected a first group of 4 massive star-forming regions 
among the northern hemisphere sources observed with the 
Effelsberg 100-m telescope where 6.7-GHz CH3OH maser 
Zeeman-splitting was measured (Vlemmings 2008). The sources 
were selected based on their peak flux density to allow 
potential detection of Zeeman-splitting in several individual 
CH3OH maser features. 

2.1. W51 

The W51 complex is one of the most luminous star-forming 
regions in the first quadrant of the Galactic plane. The high 
luminosity comes from a large number of O-type stars that 
are within the molecular cloud (Bieging 1975; Carpenter & 
Sanders 1998). Kundu & Velusamy (1967) found that W51 is 
divided in four different regions: W51A, W51B, W51C, and 
W51D. These regions are composed of several sub-regions, 
for instance W51A consists of G49.4-0.3 and G49.5-0.4 (e.g., 
Kang et al. 2010). Moreover, G49. 5-0.4 comprises the sources 
W51 -Main/South (it contains the UCHii regions W51-el and 
W51-e2; e.g., Gaume et al. 1993), W51-IRS1, and W51-IRS2 
(it contains W51-d2; e.g.. Lacy et al. 2007). 

W51-el/e2 

W51-el and W51-e2 iVx^ - 59 kms"') are the brightest 
molecular cores located at a distance of 5.41 ^qj^ kpc (Sato 
et al. 2010) in the eastern edge of W51. Their masses are 
-150 Mo and -110 Mq, respectively (Zhang & Ho 1997). 
Observations of molecular lines showed evidence for infalling, 
or accreting, gas with a possible rotation around W51-e2 (e.g., 
Zhang et al. 1998; SolHns et al. 2004; Keto & Klaassen 2008). 
Keto & Klaassen (2008) suggested a possible bipolar outflow 
along the North- West and South-East direction from W51-e2 
(PA a: 150°), which appears to be along the rotation axis of the 
ionized disc. 

OH maser polarization observations reveal that W51-e2 
has two Zeeman pairs implying a magnetic field strength 
of about -20 mG (Argon et al. 2002). Vlemmings (2008) 
measured a Zeeman-splitting of the 6.7-GHz CH3OH maser 
AVz = 0.72 + 0.04 ms"'. SubmilHmeter Array (SMA) observa- 
tions revealed that the inferred B^ morphology is hourglass-like 
near the collapsing core of W51-e2, with its pinched direction 
parallel to the direction of the ionized accretion (Tang et al. 
2009). 



W51-IRS2 

W51-1RS2 is the most luminous massive star-forming region in 
the Milky Way (Erickson & Tokunaga 1980) and it is located 
at a distance of 5.1!f^ kpc (Xu et al. 2009). W51-IRS2 hosts 
a young 05-type protostar and an UCHii region called W51- 
d2 (e.g., Gaume et al. 1993; Zapata et al. 2008). W51-d2 is lo- 
cated West w.rt. the H2O maser complex called W5 1-North at a 
projected distance of about 17000 AU. The source is associated 
with NH3 and CH3OH masers as well as H2O masers (Gaume et 
al. 1993; Minier et al. 2001; Eisner et al. 2002; Phillips & van 
Langevelde 2005). 



2.2. \N48 

The W48 complex is located in the Carina-Sagittarius spiral arm 
at a distance of 3.27+[5^^ kpc (Zhang et al. 2009). Onello et al. 
(1994) showed that W48 consists of five Hii regions: W48A- 
E. W48A hosts a cometary UCHii region called G35.2-1.74 
that is believed to be a site of massive star formation (Wood & 
Churchwell 1989; Roshi et al. 2005), and two infrared sources 
IRS 1 and IRS 2 (Zeilik & Lada 1978). CH3OH masers are 
associated with G35.2-1.74 (Caswell et al. 1995; Minier et al. 
2000). The CH3OH maser region shows a ring-like structure of 
~ 200x400 mas along a direction North-South. The mid-infrared 
source MIR 3 corresponds to G35.2-1 .74 and it is quite predomi- 
nant in the mid-infrared images of De Buizer et al. (2005). MIR 1 
coincides with the OH and H2O masers site. 

Vlemmings (2008) measured a Zeeman-splitting of the 
CH3OH masers of 0.32+0.02 ms"'. SubmilHmeter polarimetry 
observations with the Submillimetre Common User Bolometer 
Array (SCUBA) revealed two cores, the main and bright core 
in the East, W48main, and the small core in the West, W48W 
(Curran et al. 2004). The two cores appear to be connected by a 
ridge of dust and gas. The degree of polarization is lower (0.5%) 
towards the intensity peak of W48main, which is close to G35.2- 
1.74, and is higher towards the edges. Moreover, the magnetic 
field orientation is complicated and near the CH3OH maser re- 
gion, towards the centre of W48main, the magnetic field has a 
NW-SE orientation. Whereas, in W48W the magnetic field has 
an ordered orientation North-South. This changing in the ori- 
entation has been explained by twisting of the magnetic field 
towards the centre of W48main (Curran et al. 2004). 



2.3. IRAS 18556+01 38 

IRAS 18556H-0138 is at a distance of 2.19+[;-^4 kpc (Zhang et al. 
2009) and has a systemic velocity Vsys ~ 34 kms"' (Matthews 
et al. 1984). This IRAS source is associated with the mas- 
sive star-forming region G35.2-0.74N that contains an embed- 
ded BO. 5 star with a luminosity of ~ 10"* L© and a stellar 
mass of ~ 15 M^ (Dent et al. 1985; Gibb et al. 2003). A 
collimated CO-molecular bipolar outflow (PA= 58°, Gibb et 
al. 2003) was detected from G35.2-0.74N (e.g.. Dent et al. 
1985). The outflow has a size of about 2 pc and its blue- 
(25.3 kms-i< V^itf_ip,As < 30.3 kms"') and red-shifted 
(37.8 km s"' < VJu'Jf_iRAs < 43.8 km s"') parts are located South- 
West and North-East w.rt. the protostar (Lopez-Sepulcre et al. 
2009), respectively. A radio/infrared jet was detected with a po- 
sition angle PA~ 0°(Fuller et al. 2001; De Buizer 2006) sug- 
gesting that the outflow is driven by the precessing jet (e..g.; 
Marti et al. 1993). In addition, a flattened structure rotating 
perpendicularly (PA= 148°) to the outflow axis was also ob- 
served (e.g., Lopez-Sepulcre et al. 2009). Lopez-Sepulcre et al. 
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(2009) measured a velocity range of the rotating structure be- 
tween 33.45 kms-i (NW) and 34.75 kms"' (SE), and a SE-NW 
size of ~ 30", which at 2. 19 kpc correspond to about 66000 AU 
(~ 0.3 pc). 

Hutawarakorn & Cohen (1999) measured magnetic field 
strengths between -2.5 mG, North of G35.2-0.74N, and 
+5.2 mG, South by studying the Zeeman-splitting of the OH 
masers. They also found an orientation of the magnetic field par- 
allel to the outflow, even if Faraday rotation could play an im- 
portant role. Using Effelsberg observations Vlemmings (2008) 
measured a Zeeman-splitting of the CH3OH maser emission 
AVz = 0.81+ 0.04 ms-^ 

2.4. W3(0H} 

W3(OH) is the most studied UCHii region in our Galaxy. It is 
located in the Perseus spiral arm at a distance of 1.95+0.04 kpc 
(Xu et al. 2006). Infrared and radio continuum observations in- 
dicated that the central ionizing star is a 07 star with a mass 
of about 30 Mo (Scott 1981; Campbell et al. 1989). Moscadelli 
et al. (1999, 2010) derived accurate proper motion of the 12- 
GHz CH3OH masers of which the most intense concentrate 
in a small area towards the North (called the northern clump) 
of the UCHii region. A second clump is located South, and 
between the two clumps an isolated single spot maser has 
been found. This North-South filamentary structure of 12-GHz 
CH3OH masers is also shown by the 6.7-GHz CH3OH and the 
OH masers (Harvey-Smith & Cohen 2006). The brightest 12- 
GHz CH3OH masers of the northern clump are distributed along 
an axis at PA= 141° and show a regular trend of radial velocities 
with position, suggesting that they could trace a self-gravitating, 
low-mass {M x 1.5 Mq) circumstellar disk {R x 270 AU) in 
the phase of being photo-evaporated by the strong UV-radiation 
field escaping from the W3(OH) UCHii region (Moscadelli et 
al. 2010). 6.7-GHz CH3OH masers were detected for the first 
time at high angular resolution by Menten et al. (1992) and they 
coincide with the masers at 12-GHz (Moscadelli et al. 1999). In 
particular some 6.7-GHz CH3OH masers trace a similar linear 
structure as shown by the 12-GHz CH3OH masers (PA= 130°; 
Harvey-Smith & Cohen 2006). 

Polarization observations of the CH3OH masers with 
MERLIN showed that the magnetic field is parallel to the ex- 
tended filament but with a complex structure in the dominating 
northern clump (Vlemmings et al. 2006). Wright et al. (2004) 
detected linear polarization from the 1.6-GHz OH maser emis- 
sion, showing a wider spread maser polarization angles, with 
the median polarization angle (x'^a) - 104 ± 27°, that has 
been attributed to Faraday rotation along the maser path. A 
Zeeman-spHtting of the 6.7-GHz CH3OH masers AVz = 0.141 + 
0.003 m s"' was measured by Vlemmings (2008). 

3. Observations and data reduction 

The group was observed at 6.7-GHz in full polarization spec- 
tral mode with nine of the European VLBI Network' anten- 
nas (Jodrell2, Cambridge, EflFelsberg, Onsala, Medicina, Torun, 
Noto, Westerbork, and Yebes-40m), for a total observation 
time of 20 h, on November 2 (W48, IRAS 18556+0138, and 
W3(OH)) and 3 (W51), 2009 (program codes ES063A and 
ES063B). The bandwidth was 2 MHz, providing a velocity range 
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of ~ 100 kms '. The data were hardware-correlated with the 
Joint Institute for VLBI in Europe (JIVE) correlator (Schilizzi 
et al. 2001) using 1024 channels and generating all 4 polariza- 
tion combinations (RR, LL, RL, LR) with a spectral resolution 
of 1.9 kHz (~0.1 kms '). In order to detected CH3OH maser 
towards both W51-el/e2 and W51-IRS2 the data of W51 were 
correlated twice with two different correlator positions. The cor- 
relator position of W51-el/e2 was the pointing position of the 
observations. 

The data were edited and calibrated using AIPS. The band- 
pass, the delay, the phase, and the polarization calibration were 
performed on the calibrators J2202+4216 (W48 and W51) and 
3C286 (IRAS 18556+0138 and W3(OH)). Fringe-fitting and 
self-calibration were performed on the brightest maser feature 
of each star-forming region. Then the /, Q, U, RR, and LL 
cubes were imaged (2 arcsec x 2 arcsec, rms ranges from 
4 mJybeam"' to 50 mJy beam" 'depending on the source) using 
the AIPS task IMAGR. The beam-sizes were 5.5 mas x 4.5 mas 
for W51, 6.8 mas x 5.4 mas for W48, 7.5 mas x 4.8 mas 
for IRAS 18556+0318, and 8.7 mas x 4.1 mas for W3(OH). 
The Q and U cubes were combined to produce cubes of polar- 
ized intensity (Pi = VG^^^) ^nd polarization angle (x - 
1/2 X ataniUIQ)). The EVN observations were obtained close 
to VLA polarization observations^ made by the NRAO at the 
beginning of November 2009, during which the polarization an- 
gles of J2202+4216 and 3C286 were ~ -29?2 and ~ 11?5 re- 
spectively. Hence, we were able to estimate the polarization 
angles with a systemic error of no more than ~ 3°. The for- 
mal error on x are due to thermal noise. This error is given by 
a-^ = 0.5 a-pIP X 180° /;r (Wardle & Kronberg 1974), where 
P and cTp are the polarization intensity and corresponding rms 
error respectively. We were able to obtain the absolute position 
of the brightest features of W51-e2, W51-IRS2, and W3(0H) 
through fringe rate mapping using the AIPS task FRMAP. The 
absolute position errors were less than 30 mas. 



4. Analysis 

We identified the CH3OH maser features using the process de- 
scribed in Surcis et al. (2011a, 2011b). The maser features that 
showed linear polarization emission were fitted by using the 
adapted full radiative transfer method code for CH3OH masers 
described in Surcis et al. (2011a) and based on the models for 
H2O masers of Nedoluha & Watson (1992). The fit provides the 
emerging brightness temperature (TbAQ) and the intrinsic ther- 
mal linewidth (AVi). We modeled the observed linear polarized 
and total intensity maser spectra by gridding AVi from 0.5 to 
1.95 kms '. Considering TbAQ and Pi we determined the an- 
gles between the maser propagation direction and the magnetic 
field (0), if > ffcrit = 55°, where flcrit is the Van Vleck an- 
gle, the magnetic field appears to be perpendicular to the linear 
polarization vectors, otherwise it is parallel. We determined the 
Zeeman-splitting (AVz) from the cross-correlation between the 
RR and LL spectra, like was successfully used in Surcis et al. 
(2009, 201 lb) for the polarized CH3OH maser emission detected 
in W75N and NGC 7538. The dynamic range of the RR and LL 
cubes decreases close to the strongest maser emission of each 
group because of the residual calibration errors. As a result, we 
were not able to determine AVz for any CH3OH maser features 
of W3(OH) with a peak flux less than ~15 Jybeam '. 



http ://w w w. aoc . nrao. edu/~ smy ers/calibration/ 
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Fig. 1. Left panel: a view of the CH3OH maser features detected around W51-e2. Right panel: a zoom-in view of the boxed region of the left panel. 
The triangles symbols are the identified CH3OH maser features scaled logarithmically according to their peak flux density (Table 1). The maser 
LSR radial velocity is indicated by color. A 6 Jy beam"' symbol is plotted for illustration in both panels. The linear polarization vectors, scaled 
logarithmically according to polarization fraction P\, are overplotted. The two aiTows indicate the direction of the bipolar outflow (PA~ 150° ; Keto 
& Klaassen 2008). The dashed line is the best linear fit of the CH3OH maser features in the right panel (PA= 57°) 



5. Results 

In Tables 1-5 (see online material) we list all the 154 
CH3OH maser features detected towards the 4 massive star- 
forming regions observed with the EVN. The description of the 
maser distributions and the polarization results are reported for 
each source separately in Sects. 5.1-5.4. 



5.1. W51 

We detected 6.7-GHz CH3OH maser emission in the star- 
forming region W5 l-el/e2 and W5 1-IRS2. In particular, 37 fea- 
tures towards W51-e2 (named as W51E.01-W51E.37in Table 1) 
and 18 towards W51-IRS2 (named as W51N.01-W51N.18 in 
Table 2). No CH3OH maser emission (8cr = 0.1 Jy beam"') was 
detected around W5 1-e 1 . 

W51-e2 

In Fig. 1 we show the CH3OH maser features colour-coded 
according to the LSR radial velocities (50 kms"'< Visr < 
61 kms"'). The maser features near the 870 //m continuum 
peak (Tang et al. 2009) lie in an elongated structure (PA= 57°) 
perpendicular to the molecular bipolar outflow (PA~ 150°; 
Keto & Klaassen 2008) and close to the ionized accreting flow 
direction (PA~ 30°; Keto & Klaassen 2008). This suggests that 
the masers potentially trace a disc structure perpendicular to the 
bipolar outflow. However, the velocities of the CH3OH maser 
features are consistent with the velocities at the centre of the 
H53q' emission, which traces the rotating flow (Keto & Klaassen 
2008). This suggests that the masers might probe material that 
is being accreted onto the disc itself, like in Cepheus A and 
NGC7538 (Vlemmings et al. 2010; Torstensson et al. 2011; 
Surcis et al. 201 la), rather than the disc material. 

We detected linearly polarized emission (Pi - 1 .2% - 4.0%) 
in 16 CH3OH maser features that are all located close to the 
centre of W51-e2. The eiTor weighted linear polarization 
angle is (;f)w5i-e2 = 33°+16°. The full radiative transfer 
method code for CH3OH masers was able to fit 14 features. 
The intrinsic thermal linewidths and the emerging brightness 
temperatures are given in column 9 and 10 of Table 1 and 
their weighted values are (Ayi)w5i-e2 = 0-^-02 kms"' and 



(rbAQ)w5i-e2 ~ 10^ K sr, respectively. The angles between 
the maser propagation direction and the magnetic field {9) 
are almost constant for all the CH3OH maser features. The 
feature W51E.17 shows 6 = 90° indicating that the feature 
could be partially saturated. Hence, excluding this feature, 
<6')w5i-e2 = 79°+4oI. Since 6 > G^^t - 55° the magnetic field is 
most likely perpendicular to the linear polarization vectors even 
though it might also be parallel because of the relatively large 
errors in column 12 of Table 1. We measured Zeeman-splitting 
(column 1 1) in 4 CH3OH maser features. The weighted value is 
(Ayz)w5i-e2 = -1-6 ±1.2 ms"' and in absolute value is twice 
larger than that measured by Vlemmings (2008). 

W51-IRS2 

We detected 18 CH3OH maser features with peak flux density 
between 0.16 mJybeam"' and 12 mJybeam"' and with veloci- 
ties 54 kms"'< Visr < 61 kms"'. The CH3OH maser features, 
all of which are associated with the UCH 11 region d2, seem to 
lie in an arch stiTicture oriented South- West. We detected lin- 
early polarized emission only towards the brightest maser fea- 
ture (W51N.11, Pi = 1.3% and;^-w5i-iRS2 = 27°+13°). The in- 
trinsic thermal linewidth and the emerging brightness tempera- 
ture are in column 9 and 10 of Table 2 and the angle 9 (column 
13) indicates that the magnetic field is likely parallel to the linear 
polarization vector Indeed, taking into account the eiTors, 9 an- 
gles are less than G^^^i - 55°. No Zeeman-splitting was measured. 



5.2. V\/48 

In Fig. 3 all the 20 CH3OH maser features are shown, named 
in Table 3 as W48.01-W48.20. The CH3OH maser features ap- 
pear to be symmetrically distributed, with the red-shifted fea- 
tures (Visr > 43 kms"') located in the North-East and those 
blue-shifted {Vi^, < 43 kms"') in the South- West. AH the 
CH3OH maser features are associated with the submillimetre- 
core W48main detected with SCUBA by CuiTan et al. (2004) 
and are closer to MIR3 than to MIRl implying that the 
CH3OH maser features are not related to the other two maser 
species (H2O and OH). 

Linearly polarized emission was detected towards 5 
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Fig. 2. A close-up view of the CH3OH maser features detected around 
W51-IRS2. The triangles symbols are the identified CH3OH maser 
features scaled logarithmically according to their peak flux density 
(Table 2). The maser LSR radial velocity is indicated by color. A 
1 Jybeam ' symbol is plotted for illustration. The linear polarization 
vectors, scaled logarithmically according to polarization fraction Pi, are 
overplotted. 




Fig. 3. A close-up view of the CH3OH maser features detected around 
W48 (G35. 20-1.74). The triangles symbols are the identified maser 
features scaled logarithmically according to their peak flux density 
(Table 3). The maser LSR radial velocity is indicated by color. A 
30 Jybeam"' symbol is plotted for illustration. The linear polarization 
vectors, scaled logarithmically according to polarization fraction Pi, are 
overplotted. 



CH3OH maser features that show very high linear polarization 
fraction (Pi = 4.4% - 6.5%). These high values of Pi suggest 
that the features might be partially saturated. Indeed the emerg- 
ing brightness temperatures obtained by using the full radiative 
transfer method code are (for all but W48.10) greater than the 
limit below which CH3OH maser features can be considered un- 
saturated, i.e. TbAO = 2.6 x 10^ K sr (Surcis et al. 2011a). As re- 
ported in Surcis et al. (2011a), the intrinsic thermal linewidth for 
saturated masers is overestimated so the observed linewidths Avl 
might be equal to AVi (column 6 and 9 of Table 3, respectively). 
The 6 angles (column 12) are all greater than 55° so the mag- 
netic field is certainly perpendicular to the linear polarization 
vectors. Zeeman-splitting was measured only for the brightest 
feature (W48.14) and its value is AVz = 3.8 + 0.4 ms"', which 
is 10 times larger than that measured by Vlemmings (2008). 



In Table 4, named as 118556.01-118556.28, we report all the 
identified 6. 7 -GHz CH3OH maser features. In particular, they 
are associated with the source G35.20-0.74N and they can be 
divided in two groups (A and B). Group A is composed of 20 
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Fig. 4. A close-up view of the CH3OH maser features detected around 
IRAS 18556+0138 (G35.20-0.74N). The triangles symbols are the 
identified CH3OH maser features scaled logarithmically according to 
their peak flux density (Table 4). The maser LSR radial velocity is in- 
dicated by color. A 6 Jybeam"' symbol is plotted for illustration. The 
linear polarization vectors, scaled logarithmically according to polar- 
ization fraction Pi, are overplotted. The dashed line is the best linear fit 
of the CH3OH masers (PA= 140°). The two arrows indicate the direc- 
tion, but not the origin, of the coUimated CO-molecular bipolar outflow 
(PA=58°;Gibbetal. 2003). 

CH3OH maser features (27.60 km s"' < Vis, < 36.56 km s"') that 
seem to be along an arch structure (top panel of Fig. 5), while 
Group B (27 kms"'< Vis, < 29 kms"') is located about 2.3 arc- 
sec (~5000 AU at 2. 19 kpc) South-East of group A and close to 
the peak of the radio continuum emission shown in De Buizer 
(2006). The maser distribution coincides with that reported by 
Sugiyama et al. (2008). The two groups of CH3OH masers 
seem to be at the ends of a Unear structure with PA- 140°and 
a radius of about 5000 AU, which is almost perpendicular to 
the direction of the CO-outflow (PA= 58°; Gibb et al. 2003). 
This is suggested by the position angle of this linear structure 
(PA= 140°) that agrees with the position angles obtained both 
from the OH maser emission (Hutawarakom & Cohen 1999) and 
from the C0"^0(2 - 1) emission (Ldpez-Sepulcre et al. 2009). 
If we consider the velocities of the CH3OH and OH masers, of 
the large flattened structure, of the outflows, and of the source 
itself we find that the whole picture is quite different. The ve- 
locities of the two maser species are in good agreement, both 
are in the range 21 kms"'< Vmaseis < 37 kms"', and probably 
they are tracing the same gas. While the mean radial velocity 
of the large flattened structure is identical to the LSR velocity 
of the source (vG35.20 _ ^4 kms"', Hutawarakom & Cohen 
1999). The velocities of group B, which should be at the cen- 
tre of the linear structure, are about 6 kms"' blueshifted with 
respect to v'^35.20 ^j^^ jj^ j-j^g velocity range of the blue-shifted 



part of the outflow (25.3 kms < 



T/blue 

•^ outt-IRAS 



< 30.3 kms"'). 



This suggests that group B is associated with the outflow rather 
than with the linear structure. Also the velocities of group A ap- 
pear to be in agreement with the outflow velocities (25 .3 km s" ' < 
Voutf-iRAS < 43.8 kms"') although the arch-like distribution in- 
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served with MERLIN by Harvey-Smith & Cohen (2006) and 
no new group of features was observed. We do not detected 
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Fig. 5. Zoom-in view of tiie two CH3OH maser clusters of 
IRAS 18556+0138 (G35.20-0.74N). 



dicates an association with a torus structure. This torus structure 
might be related to the large flattened structure (33.45 kms"'< 
Vtorus-iRAS <34.75kms"'). 

We detected linear polarization in 5 CH3OH maser features 
(Pi = 4.8% - 11.5%) and the full radiative transfer method 
code was able to properly fit only two of them (118556.08 
and 118556.11). This is due to the likely high saturation de- 
gree of the CH3OH maser features. For the features 118556.01, 
118556.04, and 118556.27 the gridding of AVi used for fitting 
the linear polarized and total intensity spectra was insufficient 
and consequently the code provided only an upper limit for AVi. 
All the emerging brightness temperatures obtained from the fit, 
except 118556.11, reveal that the features have a high satura- 
tion degree. In particular, the reference maser feature 118556.27 
shows the highest linear polarization fraction among all the 6.7- 
GHz CH3OH maser features detected in the whole EVN group 
(Pi = 1 1.5%). In G35. 20-0.74 the magnetic field is also peipen- 
dicular to the linear polarization vectors (6* > 55°). No 6.7-GHz 
CH3OH maser Zeeman-splitting was measured in this source. 



5.4. W3(0H) 

W3(OH) is the massive star-forming region towards which we 
detected one third of all the CH3OH maser features of the First 
EVN Group. They are listed in Table 5, where they are named 
as W3OH.01-W3OH.51, and plotted in Fig. 6. They can be di- 
vided in six groups (I-VI) and their zoom-in views are shown 
in Fig. 7. Morphologically the 6.7-GHz CH3OH maser features 
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Fig. 6. A complete view of the CH3OH maser features detected around 
W3(0H). The triangles symbols are the identified CH3OH maser 
features scaled logarithmically according to their peak flux density 
(Table 5). The maser LSR radial velocity is indicated by color. A 
60 Jy beam"' symbol is plotted for illustration. The linear polarization 
vectors, scaled logarithmically according to polarization fraction Pi, are 
overplotted. 



the CH3OH maser features located West of the isolated feature 
W3OH.05 (features 22 and 23 in Harvey-Smith & Cohen 2006). 
The linear distribution of the central 6.7-GHz CH3OH maser fea- 
tures of group II (PA= 121°) is similar to those of the 12-GHz 
CH3OH masers (PA= 141°; Moscadelli et al. 2010) and of the 
OH masers (PA= 130°-140°; MoscadelU et al. 2010), which in- 
dicates that all the masers complement each other (Moscadelli 
et al. 2010). In the top-right panel of Fig. 7 the three associa- 
tions A, B, and C of 6.0-GHz OH, 12-GHz CH3OH, and 6.7- 
GHz CH3OH masers are also reported (Moscadelli et al. 2010). 
Group IV seems to lie in an arch-like structure. The northern 
clump (groups I, II, and III) shows velocities in absolute value 
greater of ~ 2-4 kms ' than the southern clump (groups IV, 
V, and VI), in particular group V has in absolute value the low- 
est velocity (Vis, ~ -42 kms"'). Along the linear structure of 
group II the distribution of the maser velocities is well ordered, 
the most blue-shifted velocities are northwestern and the most 
red-shifted are southeastern. This velocity distribution coincides 
perfectly with that of the 12-GHz CH3OH masers (Moscadelli et 
al. 2010). 

We have detected linear polarization in more than 50% of the 
CH3OH maser features (27/51) with linear polarization fraction 
ranging from 1.2% to 8.1%. The highest Pi corresponds to the 
brightest maser feature (W30H.22) that also shows the highest 
emerging brightness temperature. This value indicates that the 
feature is partially saturated like 5 other features (W3OH.05, 
W3OH.10, W30H.14, W30H.17, and W30H.18) that have 
TbAQ > 2.6 X 10^ K sr The error weighted emerging bright- 
ness temperature and the weighted intrinsic thermal linewidth of 
the unsaturated features are (rbAQ)w3(0H) ~ 2 x 10^ K sr and 
(AV'i)w3(0H) - 0.9^Q3 kms"', respectively. The error weighted 
linear polarization angle is (;k')w3(0H) - 3°+39°, while con- 
sidering only the unsaturated features (;t)w3(0H) = 15°+31°. 
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Fig. 7. Zoom-in view of the six CH3OH maser groups of W3(0H). The dashed line in panel II is the best linear fit of the central CH3OH maser 
features (PA= 121°). A, B, and C are used to label three associations of 6.0-GHz OH, 12-GHz CH3OH, and 6.7-GHz CH3OH masers. 



The values reported in column 13 of Table 5 indicate that 
for some CH3OH maser features (e.g., W30H.35) the magnetic 
field is more likely perpendicular to the linear polarization vec- 
tors, while for other features it is parallel (e.g., W3OH.09). For 
7 CH3OH maser features was possible to measure the Zeeman- 
splitting of which values are in the range -3.5 ms"'< AVz < 
3.8 ms"' (column 12 of Table 5). 



6. Discussion 

6. 1 . CHi OH maser properties 

Before discussing the magnetic fields, it is worthwhile exam- 
ining the properties of the CH3OH maser features (or simply 
CH3OH masers). As reported in Sect. 5 we observed a num- 
ber of saturated CH3OH masers in our EVN group. If we con- 
sider all the 72 CH3OH masers detected in this group and in 
NGC7538 for which was possible to determine the emerging 
brightness temperature (Surcis et al. 2011a), we can estimate the 
upper limit of linear polarization fraction for unsaturated 6.7- 
GHz CH3OH masers. In Fig 8 we report the emerging bright- 
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ness temperature as function of the linear polarization fraction. 
From this plot we find that the 6.7-GHz CH3OH masers with 
P\ < 4.5% are all unsaturated. Furthermore, from the maser the- 
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Fig. 8. The emerging brightness temperatures (Tt,AQ.) as function of the 
linear polarization fraction (Pi). The blue and red circles indicate the 
unsaturated and saturated masers, respectively, detected in NGC7538 
(Surcis et al. 20 11 a), W5 1, W48, IRAS 18556+0138, and W3(OH). The 
red arrows indicate that the TbAfl values obtained from the radiative 
transfer method code are lower limits. The red full line is the limit of 
emerging brightness temperature above which the CH3OH masers are 
considered saturated (TbAQ > 2.6 x 10' K sr; Surcis et al. 2011a), and 
the dotted line gives the lower limit to the linear polarization fraction 
for saturated masers (P\ x 4.5%). 

ory the beaming angle (AQ.) of unsaturated masers is roughly 
proportional to |ry|"', which is proportional to the brightness 
temperature T^, (Elitzur 1992). Consequently, if these masers are 
unsaturated AQ oc (Tb)"'. As shown in Surcis et al. (2011a), we 
can estimate AQ by comparing the brightness temperature T\,, 
which is given by 

■^CH,OH, (1) 



[K] 



S(y) 



[Jy] \ [mas^] / 

with TbAQ. obtained from the model. In Eq. 1 S(v) is 
the flux density, 2 the maser angular diameter, and 
fcHjOH = 13.63 ■ 10^ mas^ Jy^' K. From the Gaussian fit 
of the masers S is between 4 and 8 mas and almost all the 
CH3OH masers appear resolved. We find that AQ ranges from 
3 • 10"^ sr to 8 • 10"^ sr. In Fig. 9 aU the values of AQ as 
function of Th are reported, both in logarithmic scales. From 
a linear fit (blue line) we found that AQ oc (Tb)""'' implying 
that the CH3OH masers we have taken into account are indeed 
unsaturated. 
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Fig. 9. The beaming angle (AQ.) as function of the brightness tempera- 
tures (Ti,) for all the unsaturated CH3OH masers detected in NGC7538 
(Surcis et al. 201 la), W51, W48, IRAS 18556+0138, and W3(0H). The 
arrows indicate that the masers are unresolved, i.e. the maser angular 
size S is less than the beam-size of the observations and consequently 
AQ is smaller and T^ is larger. The blue line is the best linear fit of the 
resolved masers and the dashed black line is the theoretical proportion- 
ality An oc (T^)-^ (Elitzur 1992). 



the Lande g-factor of the 6.7-GHz CH3OH transition from which 
the Zeeman-splitting coefficient az is evaluated (Vlemmings et 
al. 201 1). If we knew the magnetic field strengths we could es- 
timate the importance of the magnetic fields by evaluating the 
ratio between thermal and magnetic energies {/3). If /3 < I, the 
magnetic field dominates the energies in high-density protostel- 
lar environment. 

Despite the uncertainty of the magnetic field strengths affects 
the jS values we can still qualitatively determine if the magnetic 
fields is important in the dynamics of the massive star-forming 
regions. Substituting AVz = az ■ By in Eq. 12 of Surcis et al. 
(2011a), where By is the magnetic field strength along the line 
of sight, and by assuming rtu^ = 10^ cm"-' and Tk ~ 200 K, we 
obtain 



/3^ 611.6 a^z^^os (6) 



lAVzl 
[ms-i] 



(2) 



Surcis et al. (20 11 a) argued that a reasonable value for az should 
be in the range 0.005 kms 'G ' < az < 0.05 kms 'G ', 
consequently from Eq. 2 we can estimate a range of /3. 

The three massive star-forming regions that host 
CH3OH masers that show circular polarization (W51, W48, and 
W30H) have j0 values between 10 ** and 10 '. This indicates 
that the magnetic fields play an important role in all the three 
regions despite we cannot exactly determine how much it is 
because of the uncertainty of Lande g-factor 

6.2.2. Magnetic field orientations 

The measurements of the linear polarization angles (x) might be 
disturbed by the rotation of the linear polarization vectors due to 
the medium between the source and the observer This rotation 
Zeeman-splitting measurements because of the uncertainty of is the foreground Faraday rotation given by Eq. 6 of Surcis et 



6.2. Magnetic fields in ttie First EVN Group 

6.2.1 . Tine importance of magnetic fields in dynamics of 
massive star formation 

We cannot determine the exact magnetic field strengths from our 
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Fig. 10. Left: modified version of Fig. 5a of Tang et al. (2009). The magnetic field (red segments) detected with the SMA (angular resolution 0".7 
that at 5.4 kpc corresponds to ~4000 AU) is superimposed on the 870 fim continuum contour map of W51-e2. The dashed black line indicates 
the direction of the ionized accreting flow by Keto & Klaassen (2008). Right: The magnetic field (red segments) detected with the SMA is 
superimposed on the 870 //m continuum contour map of W51-e2 (Tang et al. private communication). In both images the green segments, which 
are scaled logarithmically according to polarization fraction P\, mark the direction of the magnetic fields as derived from the linearly polarized 
emission of the CH3OH masers (blue triangles; angular resolution 0".001 corresponding to ~ 5 AU). 



al. (2011a). Assuming the interstellar electron density and the 
magnetic field are n^ ^ 0.012 cm"-' and By ^ 2juG (Sun et al. 
2008), Of for the 6.7-GHz CH3OH masers can be written as: 



(Df[°] = 2.26 (-^ 
\ [kpc] 



(3) 



i.e. the farther the maser source is, the larger Of is. This means 
that the foreground Faraday rotation is estimated to lie between 



(j,W5I 



12° and O 



W3(0H) 
f 



:4°. 



The quite perfect alignment of the linear polarization vectors 
of the CH3OH masers in W51-e2 suggests that the internal 
Faraday rotation is negligible. Whereas the alignments of the 
vectors in W48 and IRAS 18556+0318 additionally implies that 
the saturation state of the masers does not affect the linearly 
polarized emission. Indeed the internal Faraday rotation can 
destroy the linear polarization (Fish & Reid 2006) but we 
measured a quite high linear polarization fraction (~5%) in all 
the saturated CH3OH masers in W48 and IRAS 18556+0318. 
Hence, the contribution of the internal Faraday rotation must be 
negligible both in unsaturated and in saturated CH3OH masers 
as suggested by the alignments of the linear polarization vectors. 

W51 

Because (0)w5i-e2 = 79° > Ochi the magnetic field around 
W51-e2 is on the plane of the sky with a position angle of 
~ 120° that is perfectly perpendicular to the ionized accretion 
flow (PA~30°; Keto & Klaassen 2008) rather than to the Hnear 
distribution. The linearly polarized emission has been detected 
close to the 870 jum continuum peak and so the CH3OH masers 
are probing the magnetic field at the centre of the hourglass 
morphology mapped by Tang et al. (2009) with the SMA. In 
Fig. 10 we superimposed the results obtained by Tang et al. 
(2009) with an angular resolution 0".7 (-4000 AU) and by us 
(~ 5.4 AU). Here we can see that the magnetic field derived 
from the CH3OH masers (green segments) is consistent with the 
hourglass morphology (red segments). The negative value of the 
Zeeman-splitting indicates that the magnetic field is oriented 



towards the observer. 

In the case of W51-IRS2 ^5i_e2 = 61:*:5J° indicating 
that the magnetic field is more likely parallel to the linear 
polarization vectors than perpendicular. Consequently we have 
also here a perfect agreement with the magnetic field detected 
with the SMA at 870 //m by Tang et al. (private communica- 
tion; Fig. 10), implying a well-ordered magnetic field around d2. 

W48 

A comparison between the magnetic fields derived from the 
6.7-GHz CH3OH maser observations and from the SCUBA 
850 yum observations (Curran et al. 2004) shows that, like in 
W51, the magnetic fields are consistent towards the core at 
different angular resolution. The positive Zeeman-splitting 
(|Ayzlw48 = 3.8 ms"') indicates that the magnetic field is 
oriented away from the observer. 

IRAS 18556+0318 

Four masers of group A show linearly polarized emission from 
which the magnetic field is estimated to be oriented East- West 
(R4 = 90°). However, because of the velocity difference among 
the masers it is worthwhile to determine the orientation of 
the magnetic field separately by considering the masers two 
by two. Considering 118556.04 and 118556.11, which have 
velocities 7 kms"' greater than 118556.01 and 118556.08, the 
magnetic field is estimated to have PA=112°, while from the 
other two masers we estimated PA=83°. From group B we 
get PA=80°. In all cases the magnetic field derived from the 
6.7-GHz CH3OH masers is aligned neither with the outflow 
nor with the linear structure, but is almost perpendicular to the 
radio/infrared jet. 



6.2.3. W3(0H) 

Since we observed at higher angular resolution than Vlemmings 
et al. (2006) we are able to better determine the magnetic field 
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morphology in group II. Comparing Fig.2 of Vlemmings et al. 
(2006) and our Fig. 6 we see that the linear polarization vectors of 
the 6.7-GHz CH3OH masers are not as similar through the whole 
filamentary structure as observed with the MERLIN, probably 
because of the 90°-flip phenomenon. However, they appear well 
ordered within each group (Fig. 7) and this again indicate that 
the saturation state of the masers does not influence the linearly 
polarized emission. 

Groups III, V, and VI show error weighted linear polariza- 
tion angles (j() close to 0°. A comparison with the results ob- 
tained by Vlemmings et al. (2006), who did not detect any lin- 
early polarized emission from group V, indicates that the ori- 
entation of the linear polarization vectors in these groups are 
likely influenced by the 90° -flip phenomenon, in particular for 
group III. Althought the 90°-flip seems to be confirmed from 
the 9 angles measured by us for all the three groups (column 
14 of Table 5), in groups V and VI the magnetic field is more 
likely perpendicular to the vectors, i.e. with orientation West- 
East. Whereas in group III the magnetic field has an orienta- 
tion North-South. In group IV we see that all the features but 
W30H.18 and W30H.25 show a linear polarization angle sim- 



100 


1 1 1 


W3(0H) - 
C 






>^ 







;^ ^ 




B ^ 


100 


- 


l-^ A 






i 


200 


- 


1 - 




- A60 Jy 


bcam-^ magnetic field 




", , 1 


, , , 1 , , , , 1 , , , , 1 , ,- 



200 100 -100 

Right Ascension offset (mas) 

Fig. 11. Modified version of panel II of Fig. 7. The green segments, 
which are scaled logarithmically according to polarization fraction Pi, 
mark the direction of the magnetic field as derived from the linearly 
polarized emission of the CH3OH masers in the northern clump of 
W3(0H). 



ilar to that measured by Vlemmings et al. (2006). The different 
angles of W30H.18 and W30H.25 could also be due to the 90°- 
flip and consequently we find that here, like in group III, the 
magnetic field is oriented almost North-South. This orientation 
of the magnetic field is also determined from the isolated feature 
W3OH.05 of which linear polarization angle agrees quite per- 
fectly with that measured with MERLIN. While the magnetic 
field of groups III and VI has an orientation along the filamen- 
tary structure, the magnetic field morphology in group II is com- 
pletely different. 

For all the CH3OH masers of the group II-B that have x 
close to 0° (W30H.22, W30H.23, W30H.35, W30H.37, and 
W30H.39) we find, considering the errors, that their 6 angles 
are likely greater than 55°, indicating that the magnetic field 
is perpendicular to the large-scale filamentary structure and al- 
most aligned with the linear distribution of the B association. 
This orientation of the magnetic field is also suggested by the 



masers that have linear polarization vectors parallel to the linear 
structure (W30H.24 and W30H.28) that show 9^in < 55°, i.e. 
magnetic field is likely parallel to them. As show in Fig. 11, the 
CH3OH masers of A and C associations (W30H.27, W30H.31, 
W30H.41, W30H.42, and W30H.43) seems to trace a magnetic 
field that is changing its morphology from a North-South orien- 
tation to a West-East close to the linear structure (B association). 
From the Zeeman-splitting measurements we can determine 
the pointing direction of the magnetic fields. Let consider first 
group II. Since we measured both positive (W30H.22) and neg- 
ative (W30H.35 and W30H.37) Zeeman-splittings at the North- 
West end and at the South-East end of the linear distribution 
respectively, the magnetic field must have a counterclockwise 
direction. This direction is opposite to the rotation of the lin- 
ear distribution, which is clockwise as derived from the OH and 
CH3OH masers velocities. If the masers along the linear distri- 
bution trace a disc that is rotating clockwise (M » 1 .5 M0 and 
R ^ 270 AU; Moscadelli et al. 2010), we have in W3(OH) a 
scenario similar to that of NGC7538 (Surcis et al. 2011a). The 
magnetic field is on the disc (in NGC7538 on the torus) rotat- 
ing in opposite direction to the disc. The Zeeman-splitting mea- 
surements of the other groups suggest that the magnetic field is 
twisted along the filamentary structure. 

7. Conclusions 

We observed a group of 4 massive star-forming regions at 
6.7-GHz in full polarization spectral mode with the EVN 
to detect linear and circular polarization emission from 
die CH3OH masers. We detected a total of 154 6.7-GHz 
CH3OH masers one third of which only towards W3(OH). 
We detected linearly polarized emission in all the regions 
of the EVN group, and in 3 of them we also measured 
Zeeman-splitting. From our observations we have found that 
the magnetic field is well ordered in all the massive star- 
forming regions. In W51-e2 the magnetic field derived from the 
CH3OH masers is at the centre and along the axis of symmetry 
of the hourglass morphology observed with the SMA at 870 //m. 
In W3(OH) we determined two different magnetic fields, one 
twisted along the filamentary structure and the other one in the 
northern clump on the likely low-mass rotating disc probed by 
the linear distribution of the CH3OH and OH masers with a 
counterclockwise orientation. In the case of IRAS 18556+03 18 
(G35.20-0.74N) the magnetic field seems to be perpendicular 
to the radio/infrared jet rather than along the outflow or along 
a possible torus structure. The magnetic fields are in perfect 
agreement with their morphology at lower angular resolution 
determined with the SMA and SCUBA (W51-e2, W51-d2, 
and W48) and they play an important role in all the massive 
star- forming regions of the First EVN Group. 

In addition, we have also estimated a limit of linear polariza- 
tion fraction of 4.5% below which the 6.7-GHz CH3OH masers 
are unsaturated. 
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Table 1. All 6.7-GHz methanol maser features detected in W51~e2. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


Maser 


RA" 
offset 


Dec" 
offset 


Peak flux 
Density(I) 


Vkr 


AVL 


P\ 


X 


hV," 


T'bAn'' 


AVz'- 


(f 




(mas) 


(mas) 


(Jy/beam) 


(km/s) 


(km/s) 


(%) 


(°) 


(km/s) 


(log K sr) 


(m/s) 


(°) 


W51E.01 


^791.826 


-357.414 


0.78 ± 0.02 


51.54 


0.24 


- 


- 


- 


- 


- 


- 


W51E.02 


-741.361 


313.950 


0.24 ±0.01 


53.03 


0.22 


- 


- 


- 


- 


- 


- 


W51E.03 


-396.633 


-535.208 


2.31 ±0.02 


51.45 


0.42 


- 


- 


- 


- 


- 


- 


W51E.04 


-368.770 


-501.720 


0.47 ± 0.01 


52.42 


0.22 


- 


- 


- 


- 


- 


- 


W51E.05 


-202.860 


-304.199 


1.11 ±0.04 


51.10 


0.33 


- 


- 


- 


- 


- 


- 


W51E.06 


-198.317 


-300.964 


9.73 ± 0.04 


51.80 


0.27 


1.6 ±0.2 


40 ±12 


l-3!Si 


8-9!Si 


- 


79!i|! 


W51E.07 


-193.719 


-310.570 


5.45 ± 0.03 


52.15 


0.27 


2.2 ±0.3 


35 ±49 


1 9+o:i 

-0.2 


9.0!»-^ 


- 


80!^° 


W51E.08 


-185.520 


-311.840 


0.13 ±0.01 


52.94 


0.33 


- 


- 




- 


- 


- 


W51E.09 


-116.608 


558.552 


3.84 ±0.16 


59.35 


0.21 


- 


- 


- 


- 


- 


- 


W51E.10 


-20.773 


16.247 


2.05 ± 0.06 


60.14 


0.17 


- 


- 


- 


- 


- 


- 


W51E.11 


-16.176 


160.938 


2.02 ± 0.05 


59.88 


0.25 


- 


- 


- 


- 


- 


- 


W51E.12 


-7.368 


-65.472 


22.72 ± 0.03 


58.30 


0.27 


2.6 ±0.6 


20 ±6 


0.8!JJ] 


9 [+0.5 

8.9+1^ 


- 


79!|> 


W51E.13 


-2.936 


-65.970 


65.78 ± 0.20 


57.86 


0.36 


1.8 ±0.4 


21 ±5 


1.7+»:' 


-0.4 ±0.1 


77+!] 


W51E.14 
W51E.15 




0.443 



12.135 


217.32 ±0.23 
26.63 ± 0.07 


59.26 
59.00 


0.33 
0.25 


1.4 ±0.8 
1.4 ±0.4 


46 ±9 
29 ±18 


-0.] 




-0.7 ±0.1 


67^ 
74!i? 


W51E.16 


0.831 


-76.960 


2.28 ± 0.03 


58.39 


0.27 


- 


- 






- 




W51E.17 


1.219 


-14.587 


50.54 ± 0.05 


59.97 


0.38 


1.9 ±0.4 


35 ±15 


^■^"-il 


^Kl 


- 


90!|] 


W51E.18 


4.819 


9.281 


27.56 ±0.16 


59.35 


0.28 


1.2 ±0.2 


-2 ±32 






-2.4 ±0.1 




W51E.19 


10.581 


38.508 


4.09 ± 0.05 


59.97 


0.27 


1.5 ± 1.2 


52 ±17 


o-7!S:l 


8.8!?;J 


- 


68!i? 


W51E.20 


11.412 


-24.912 


2.94 ± 0.05 


60.05 


0.25 


- 


- 






- 




W51E.21 


27.975 


-43.551 


3.15 ±0.01 


60.32 


0.23 


2.3 ± 1.6 


39 ±46 


- 


- 


- 


- 


W51E.22 


56.393 


-143.511 


2.09 ± 0.04 


58.91 


0.19 


- 


- 


- 


- 


- 


- 


W51E.23 


67.250 


-145.042 


15.37 ±0.05 


58.74 


0.28 


2.0 ±0.5 


29 ±11 


i-i;Si 


8-9!2? 


- 


78-3? 


W51E.24 


69.300 


-158.505 


9.48 ± 0.07 


59.00 


0.30 


4.0 ± 1.6 


31±6 


0.6!Si 


9.4+»| 


- 


84+5^ 


W51E.25 


70.851 


-154.001 


8.19 ±0.03 


58.82 


0.26 


2.5± 1.1 


26 ±8 


o-7!s;i 


9.1!?;] 


- 


76!iS 


W51E.26 


168.514 


-12.121 


4.27 ± 0.03 


58.56 


0.51 


- 


- 






- 




W51E.27 


196.655 


99.159 


1.82 ±0.05 


59.80 


0.23 


- 


- 


- 


- 


- 


- 


W51E.28 


204.632 


113.516 


5.67 ± 0.02 


60.58 


0.34 


2.9 ±0.6 


19±8 


- 


- 


- 


- 


W51E.29 


262.520 


43.556 


2.26 ± 0.07 


58.91 


0.29 


- 


- 


- 


- 


- 


- 


W5IE.30 


268.004 


38.809 


1.99 ±0.03 


58.82 


0.27 


- 


- 


- 


- 


- 


- 


W51E.31 


271.716 


40.205 


3.30 ± 0.03 


58.47 


0.40 


3.5± 1.1 


63 ±5 


i.4!Si 


S-S^o^ 


- 


84!«„ 


W51E.32 


276.535 


31.513 


8.11 ±0.06 


57.86 


0.48 


1.6 ±0.5 


32 ±23 


0-<l 


-2.7 ±0.1 


75!J? 


W51E.33 


616.831 


909.365 


0.57 ± 0.07 


51.01 


0.25 


- 


- 






- 




W51E.34 


620.709 


902.697 


0.66 ± 0.02 


50.04 


0.40 


- 


- 


- 


- 


- 


- 


W51E.35 


661.979 


343.477 


0.62 ± 0.03 


58.47 


0.22 


- 


- 


- 


- 


- 


- 


W51E.36 


909.542 


98.515 


0.69 ± 0.03 


58.82 


0.28 


- 


- 


- 


- 


- 


- 


W51E.37 


923.834 


88.762 


1.29 ±0.04 


58.47 


0.30 


- 


- 


- 


- 


- 


- 



Notes. <"> The reference position is 0-2000 = 19''23"43'.942 and 6200Q - +14°30'34'.'550. <*' The best-fitting results obtained by using a model 
based on the radiative transfer theory of methanol masers for F -1- F^ = 1 (Vlemmings et al. 2010, Surcis et al. 201 lb). The errors were determined 
by analyzing the full probability distribution function. *^' The Zeeman- splittings are determined from the cross-correlation between the RR and 
LL spectra. ** The angle between the magnetic field and the maser propagation direction is determined by using the observed P\ and the fitted 
emerging brightness temperature. The errors were determined by analyzing the full probability distribution function. 
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Table 2. All 6.7-GHz methanol maser features detected in W51-IRS2. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Maser 


RA" 
offset 


Dec- 
offset 


Peak flux 
Density(I) 


Vi,, 


AVL 


P\ 


X 


AVi* 


nMt 


AVz'' 


ff 




(mas) 


(mas) 


(Jy/beam) 


(km/s) 


(km/s) 


(%) 


o 


(km/s) 


(log K sr) 


(m/s) 


O 



W51N.01 


^85.860 


-28.289 


0.27 ± 0.01 


54.52 


0.31 


W51N.02 


-83.367 


-22.108 


0.20 ± 0.01 


54.96 


0.28 


W51N.03 


-53.067 


-11.479 


0.33 ± 0.04 


57.94 


0.28 


W51N.04 


-45.312 


8.416 


0.90 ± 0.02 


57.51 


0.25 


W51N.05 


-33.569 


-351.122 


0.61 ±0.02 


59.70 


0.18 


W51N.06 


-25.536 


-312.023 


0.47 ± 0.03 


59.88 


0.28 


W51N.07 


-12.962 


13.928 


0.74 ± 0.02 


55.84 


0.25 


W51N.08 


-9.140 


-108.629 


1.02 ±0.03 


56.28 


0.29 


W51N.09 


-6.813 


-2.965 


1.93 ±0.01 


56.45 


0.21 


W51N.10 


-5.152 


-16.856 


1.12 ±0.03 


56.28 


0.35 


W51N.11 








12.03 ±0.07 


56.19 


0.35 


W51N.12 


1.496 


18.406 


0.95 ± 0.01 


56.63 


0.36 


W51N.13 


5.096 


-6.713 


1.25 ±0.01 


55.75 


0.30 


W51N.14 


5.373 


-161.138 


2.42 ± 0.05 


59.79 


0.24 


W51N.15 


10.192 


6.065 


0.57 ± 0.01 


56.98 


0.31 


W51N.16 


12.630 


-129.178 


0.48 ± 0.02 


60.84 


0.26 


W51N.17 


24.761 


-149.938 


0.16 ±0.01 


60.23 


0.26 


W51N.18 


64.644 


26.281 


0.20 ± 0.02 


55.75 


0.28 



1.3 ±1.1 27 ±13 1.8+i!-| 8.7+?? - 61+;! 



Notes. *"' The reference position is 0-2000 = 19''23™39^.818 and ^2000 = +14°3r5"045. **' The best-fitting results obtained by using a model based 
on the radiative transfer theory of methanol masers for F -1- F^ = 1 (Vlemmings et al. 2010, Surcis et al. 2011b). The errors were determined 
by analyzing the full probability distribution function. *''' The Zeeman- splittings are determined from the cross-correlation between the RR and 
LL spectra. *''' The angle between the magnetic field and the maser propagation direction is determined by using the observed P\ and the fitted 
emerging brightness temperature. The errors were determined by analyzing the full probability distribution function. 



Table 3. All 6.7 GHz methanol maser features detected in W48. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Maser 


RA 
offset 


Dec 
offset 


Peak flux 
Density(I) 


Vk, 


AVL 


P\ 


X 


AV," 


T-bAQ" 


AVz" 


8' 




(mas) 


(mas) 


(Jy/beam) 


(km/s) 


(km/s) 


(%) 


(°) 


(km/s) 


(log K sr) 


(m/s) 


(°) 



W48.01 


-138.785 


-87.692 


19.00 ± 0.40 


43.78 


0.22 


- 


- 


W48.02 


-110.296 


-57.672 


6.84 ± 0.06 


44.13 


0.19 


- 


- 


W48.03 


-109.953 


-128.784 


14.07 ± 0.26 


41.41 


0.23 


- 


- 


W48.04'' 


-101.143 


-70.511 


60.05 ± 1.46 


42.91 


0.29 


5.1 ±0.2 


23 ±8 


W48.05 


-99.026 


-134.666 


23.28 ±0.95 


40.89 


0.29 


- 


- 


W48.06 


-94.621 


-117.296 


8.38 ±0.14 


42.03 


0.24 


- 


- 


W48.07 


-90.502 


-158.745 


78.93 ± 1.40 


41.50 


0.22 


- 


- 


W48.08 


-89.072 


-140.818 


1.58 ±0.02 


40.18 


0.35 


- 


- 


W48.09 


-82.093 


-128.763 


3.82 ±0.09 


40.18 


0.35 


- 


- 


W48.10 


-81.635 


-122.286 


202.87 ± 6.56 


42.47 


0.26 


4.4 ±1.3 


22 ±6 


W48.ll 


-81.235 


-157.660 


47.04 ± 0.08 


40.97 


0.31 


- 


- 


W48.12 


-77.974 


-200.914 


22.79 ±0.18 


41.76 


0.37 


- 


- 


W48.13 


-75.171 


-192.154 


11.58 ±0.26 


41.41 


0.17 


- 


- 


W48.14'' 








294.68 ±1.11 


44.49 


0.24 


4.6 ±1.5 


12 ±6 


W48.15 


13.158 


-65.594 


2.48 ±0.12 


43.34 


0.15 


- 


- 


W48.16 


27.631 


-0.439 


26.15 ±0.41 


44.57 


0.26 


- 


- 


W48.17 


27.631 


-29.171 


22.64±1.11 


44.49 


0.12 


- 


- 


W48.18 


34.267 


-20.596 


3.05 ±0.17 


45.36 


0.24 


- 


- 


W48.19'' 


47.082 


-62.460 


101.52 ± 1.30 


43.70 


0.25 


5.4 ± 0.6 


16 ±4 


W48.20'' 


58.352 


6.796 


32.74 ±0.23 


44.84 


0.25 


6.5 ± 0.6 


27 ±1 



7+"-^ 9 4+"-^ - 84+' 



1.0+"' 9.5+?? 3.8 ±0.4 73+" 



0.7+?? 9.6+?;! - 84+* 



9 7+"--' - 86+^ 



Notes. *"' The best-fitting results obtained by using a model based on the radiative transfer theory of methanol masers for F -1- F^ = 1 (Vlemmings 
et al. 2010, Surcis et al. 2011b). The errors were determined by analyzing the full probability distribution function. *** The Zeeman-splittings 
are determined from the cross-correlation between the RR and LL spectra. *^* The angle between the magnetic field and the maser propagation 
direction is determined by using the observed P] and the fitted emerging brightness temperature. The errors were determined by analyzing the 
full probability distribution function. '^ Because of the degree of the saturation of these CH3OH masers TbAfi is underestimated, AVj and 6 are 
overestimated. 



Surcis et al.: EVN group: magnetic field., Online Material p 3 
Table 4. All 6.7 GHz methanol maser features detected in IRAS18556+01 (G35.20~0.74N). 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Maser 


Group 


RA 

offset 


Dec 

offset 


Peak flux 
Density(I) 


Vi„ 


AVL 


Pi 


X 


AV," 


TbAn" 


AVz* 


9' 






(mas) 


(mas) 


(Jy/beam) 


(km/s) 


(km/s) 


(%) 


n 


(km/s) 


(log K sr) 


(m/s) 


(°) 


118556.01" 


A 


-1587.159 


1817.554 


14.71 ±0.34 


29.89 


0.23 


8.6 ±1.3 


-9 ±3 


<0.5 


10.0«;^ 


- 


w\t 


118556.02 


A 


-1563.866 


1764.340 


1.05 ±0.10 


28.57 


0.42 


- 


- 


- 




- 




118556.03 


A 


-1561.521 


1881.611 


0.85 ± 0.03 


33.22 


0.30 


- 


- 


- 


- 


- 


- 


118556.04" 


A 


-1558.818 


1630.295 


4.66 ± 0.03 


36.56 


0.26 


9.1 ± 1.4 


22 ±28 


<0.5 


9.^11 


- 


90-11 


118556.05 


A 


-1559.276 


1826.489 


4.47 ± 0.03 


27.95 


0.24 


- 


- 


- 




- 




118556.06 


A 


-1553.628 


1790.443 


0.32 ± 0.03 


27.95 


0.26 


- 


- 


- 


- 


- 


- 


118556.07 


A 


-1551.154 


1729.794 


0.69 ± 0.02 


27.60 


0.28 


- 


- 


- 


- 


- 


- 


118556.08" 


A 


-1550.882 


1862.656 


12.13 ±0.08 


29.36 


0.29 


6.5 ± 0.7 


2±3 


l-0-a5 


9 7+0.4 

-0.2 


- 


84!t, 


118556.09 


A 


-1549.567 


1674.839 


1.31 ±0.02 


30.41 


0.30 


- 


- 






- 




118556.10 


A 


-1548.166 


1827.927 


0.52 ±0.02 


27.25 


0.22 


- 


- 


- 


- 


- 


- 


118556.11 


A 


-1548.080 


1604.614 


12.81 ±0.11 


35.42 


0.30 


4.8 ±1.1 


22 ±18 


o-7!l!:5 


94+0.4 


- 


M% 


118556.12 


A 


-1538.571 


1715.225 


0.14 ±0.01 


33.49 


0.30 


- 


- 






- 




118556.13 


A 


-1535.311 


1903.782 


0.97 ± 0.03 


33.22 


0.20 


- 


- 


- 


- 


- 


- 


118556.14 


A 


-1527.890 


1870.255 


0.56 ±0.02 


33.40 


0.14 


- 


- 


- 


- 


- 


- 


118556.15 


A 


-1519.568 


1904.572 


1.74 ±0.01 


26.90 


0.43 


- 


- 


- 


- 


- 


- 


118556.16 


A 


-1513.805 


1912.239 


1.36 ±0.01 


26.20 


0.31 


- 


- 


- 


- 


- 


- 


118556.17 


A 


-1501.437 


1774.460 


0.22 ±0.01 


30.94 


0.21 


- 


- 


- 


- 


- 


- 


118556.18 


A 


-1466.290 


1754.551 


0.20 ±0.01 


32.34 


0.37 


- 


- 


- 


- 


- 


- 


118556.19 


A 


-1454.021 


1765.602 


0.14 ±0.01 


31.73 


0.27 


- 


_ 


- 


- 


- 


- 


118556.20 


A 


-1293.200 


1994.118 


4.37 ±0.19 


29.10 


0.21 


- 


- 


- 


- 


- 


- 


118556.21 


B 


-14.056 


20.645 


1.41 ±0.02 


27.95 


0.40 


~ 


- 


- 


- 


- 


- 


118556.22 


B 


-14.013 


28.847 


3.86 ±0.05 


28.31 


0.23 


_ 


- 


- 


- 


- 


- 


118556.23 


B 


-12.669 


72.609 


1.45 ±0.05 


28.31 


0.65 


- 


- 


- 


- 


- 


- 


118556.24 


B 


-11.039 


17.769 


0.39 ±0.01 


27.43 


0.25 


- 


- 


- 


- 


- 


- 


118556.25 


B 


-5.863 


-4.211 


0.31 ±0.01 


27.34 


0.19 


- 


- 


- 


- 


- 


- 


118556.26 


B 


-3.117 


-277.435 


1.58 ±0.05 


28.39 


0.21 


- 


- 


- 


- 


- 


- 


118556.27" 


B 








104.81 ±1.05 


28.57 


0.26 


11.5 ±1.9 


-10±5 


<0.5 


10-8-a7 


- 


85^5 


118556.28 


B 


49.389 


48.145 


0.89 ± 0.03 


27.69 


0.18 


- 


- 


- 




- 


- 



Notes. *"' The best-fitting results obtained by using a model based on the radiative transfer theory of methanol masers for F -f F^ = 1 (Vlemmings 
et al. 2010, Surcis et al. 2011b). The errors were determined by analyzing the full probability distribution function. *** The Zeeman-splittings 
are determined from the cross-correlation between the RR and LL spectra. *^* The angle between the magnetic field and the maser propagation 
direction is determined by using the observed P\ and the fitted emerging brightness temperature. The errors were determined by analyzing the 
full probability distribution function. *''' Because of the degree of the saturation of these CH3OH masers TbAfi is underestimated, AVi and 6 are 
overestimated. 
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Table 5. All 6.7-GHz methanol maser features detected in W3(OH). 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Maser 


Group 


RA" 
offset 


Dec" 

offset 


Peak flux 
Density(I) 


Vw 


Avl 


P\ 


X 


hVf 


T-bAn* 


AVz'' 


e^ 






(mas) 


(mas) 


(Jy/beam) 


(km/s) 


(km/s) 


(%) 


O 


(km/s) 


(log K sr) 


(m/s) 


(°) 


W3OH.01 


III 


-834.479 


148.296 


11.58 ±0.02 


-46.55 


0.20 


1.2 ±0.4 


-2 ±10 


0.9^1 


S-S-a's 


- 


71!'!; 


W3OH.02 


III 


-817.800 


153.654 


2.26 ±0.01 


-47.07 


0.18 


- 


- 






- 




W3OH.03 


III 


-802.674 


79.334 


27.78 ±0.28 


-45.76 


0.22 


3.2 ±0.8 


12 ±69 


l-l-a2 


9.2;g;« 


- 


69-42 


W3OH.04 


VI 


-383.492 


-1764.656 


7.78 ± 0.04 


-42.24 


0.19 


- 


- 






- 


- 


W3OH.05'- 


- 


-162.497 


-691.895 


33.61 ±0.55 


-44.00 


0.18 


5.6 ± 1.7 


-69 ±5 


0-7tl 


9.6!?;^ 


- 


83!l4 


W3OH.06 


VI 


-122.187 


-1749.006 


27.32 ± 0.05 


-42.33 


0.21 


- 


- 






- 




W3OH.07 


VI 


-153.206 


-1767.649 


54.81 ±0.34 


-42.60 


0.20 


1.9 ±0.6 


35 ±14 


0.9t| 


9.0:g;5 


- 


79+18 
-34 


W3OH.08 


VI 


-149.964 


-1753.939 


4.82 ± 0.04 


-42.24 


0.19 


- 


- 






- 




W3OH.09 


VI 


-132.122 


-1771.915 


235.45 ± 0.33 


-42.60 


0.19 


2.5 ± 1.2 


15 ±33 


0.9t| 


9.2!g-^ 


- 


6TJ,1 


W3OH.10' 


VI 


-131.591 


-1742.706 


60.79 ± 0.34 


-42.60 


0.18 


5.3 ± 2.7 


-Il±7 


0-8tl 


9.6!?;^ 


- 


1T}1 


W30H.I1 


VI 


-123.541 


-1772.516 


212.03 ±0.33 


-42.60 


0.22 


- 


- 






-3.5 ±0.2 




W30H.I2 


VI 


-114.142 


-1773.373 


1 11.69 ±0.37 


-42.95 


0.20 


- 


- 


- 


- 


- 


- 


W30H.13 


VI 


-111.926 


-1747.892 


22.93 ± 0.03 


-42.07 


0.20 


3.1 ±0.3 


5 ±38 


o.sti 


9.2tJ 


- 


81+?4 


W30H.14'- 


IV 


-59.565 


-1173.670 


182.60 ±0.64 


-43.65 


0.26 


5.9 ± 1.9 


-41 ±44 


l-l^az 


9.6!?;' 


- 


75!^5 


W30H.15 


II 


-57.240 


-157.034 


20.812 ±0.59 


-44.79 


0.33 


- 


- 






- 




W30H.16 


IV 


-46.725 


-1160.805 


233.50 ± 0.38 


-42.95 


0.32 


3.9 ±0.5 


-39 ± 86 


1.5t5 


9.3t5 


- 


74!JJ 


W30H.17'' 


IV 


-42.741 


-1165.085 


1 10.50 ±0.29 


-43.12 


0.25 


5.6 ±3.9 


-75 ± 41 




9.6!fJ 


- 


75J 


W30H.18'' 


IV 


-40.408 


-1160.982 


65.18 ±0.49 


-43.12 


0.29 


5.0 ±2.6 


2 ±85 


1 2+''-' 

-0.2 


9.5:?;^ 


- 


74!i^ 


W30H.19 


II 


-35.097 


-11.770 


69.09 ± 0.62 


-44.79 


0.25 


- 


- 






- 




W3OH.20 


II 


-30.01 1 


-17.119 


46.52 ± 0.72 


-44.70 


0.25 


- 


- 


- 


- 


- 


- 


W30H.21 


IV 


-29.630 


-1310.143 


57.91 ±0.52 


-43.91 


0.31 


1.5± 1.1 


62 ±49 


1 5+°' 
LO^j 


8-8-oI 


- 


73+i8 


W30H.22'' 


II 








2051.30 ±1.42 


-45.41 


0.31 


8.1 ±0.7 


4 ±34 


I0.0!|]-5 


1.9 ±0.1 


W30H.23 


II 


11.410 


-36.799 


1099.80 ±0.81 


-43.39 


0.29 


4.1 ± 1.2 


2±41 


1-2-02 


9 4+0-6 

"•4!^! 


- 


75+" 


W30H.24 


II 


14.956 


-21.265 


913.15 ±0.92 


-44.53 


0.37 


4.4 ± 2.0 


-38 ±41 


I 7(1+0.03 


- 


W30H.25 


IV 


15.969 


-1376.160 


156.18 ±0.65 


-43.74 


0.21 


1.5 ±0.4 


7 ±23 


1 I+° ' 
-0.1 


8-8-aI 


3.8 ± 0.5 


68!^, 


W30H.26 


II 


18.653 


-39.053 


546.64 ± 0.63 


-43.65 


0.21 


- 


- 






- 




W30H.27 


II 


22.492 


-84.309 


206.85 ± 0.36 


-43.03 


0.18 


1.8 ±0.2 


63 ±35 


0.9t| 


9.0!g-5 


- 




W30H.28 


II 


24.383 


-22.999 


640.49 ± 0.70 


-44.70 


0.29 


4.1 ± 1.6 


-47 ± 19 


^■<2 


9 4+o:8 

-1.0 


- 


W30H.29 


II 


30.256 


-34.817 


239.84 ± 0.56 


-44.09 


0.27 


- 


- 






- 




W3OH.30 


II 


33.002 


-45.343 


20.09 ± 0.63 


-43.56 


0.26 


- 


- 


- 


- 


- 


- 


W30H.31 


II 


39.591 


-125.509 


74.16 ±0.50 


-43.30 


0.19 


1.2 ±0.9 


64 ±29 


0.9^1 


8.8!;-^ 


- 


65^2 


W30H.32 


II 


47.238 


-59.534 


155.18 ±0.36 


-43.03 


0.18 


- 


- 






- 




W30H.33 


II 


51.173 


-159.398 


146.66 ± 0.36 


-43.03 


0.26 


4.1 ±2.7 


6 ±42 


l-2t^ 


9A*_1l 


- 


70!i^ 


W30H.34 


II 


53.063 


-67.207 


140.91 ±0.98 


-44.53 


0.40 


- 


- 






- 




W30H.35 


II 


58.220 


-63.723 


347.43 ± 0.37 


-42.86 


0.19 


3.2± 1.1 


I ±37 


0-m 


9.2!g;^ 


-I.5±0.1 


7I!5« 


W30H.36 


II 


62.103 


-65.691 


180.85 ±0.70 


-44.35 


0.37 


- 


- 






- 




W30H.37 


II 


100.540 


-93.575 


110.11 ±0.28 


-42.51 


0.31 


3.2 ±0.8 


19 ±53 


i-4!l!i 


9-2!l!i 


-1.7 ±0.4 


76!^^ 


W30H.38 


II 


I0I.04I 


36.774 


7.17 ±0.02 


-46.37 


0.25 


- 


- 






- 




W30H.39 


II 


104.515 


-84.400 


28.12 ±0.05 


-42.33 


0.27 


2.7 ± 1.0 


-I4±6 


0-8-a? 


9-3!?;t 


- 


%\% 


W3OH.40 


I 


105.043 


540.987 


178.62 ± 0.45 


-45.14 


0.22 


- 


- 






-3.4 ±0.4 




W30H.41 


II 


106.053 


32.572 


19.21 ±0.28 


-45.76 


0.31 


2.1 ±0.3 


63 ±80 


^■'-0.2 


9-0-02 


- 


69!$^ 


W30H.42 


II 


106.329 


41.679 


2.08 ±0.01 


-46.99 


0.44 


4.4 ±0.5 


68 ±5 






- 




W30H.43 


II 


109.883 


47.815 


254.77 ± 0.69 


-44.35 


0.24 


2.6 ± 0.7 


16 ±27 


l-2tl 


"•It? 


- 


73!il 


W30H.44 


I 


117.485 


535.473 


73.74 ± 0.44 


-45.05 


0.17 


- 


- 






- 




W30H.45 


I 


170.642 


437.187 


0.72 ±0.01 


-47.07 


0.18 


- 


- 


- 


- 


- 


- 


W30H.46 


V 


399.814 


-1735.970 


4.32 ±0.11 


-42.42 


0.43 


- 


- 


- 


- 


- 


- 


W30H.47 


V 


423.769 


-1737.766 


1.25 ±0.02 


-41.89 


0.28 


- 


- 


- 


- 


- 


- 


W30H.48 


V 


437.887 


-1731.754 


19.97 ± 0.02 


-41.81 


0.78 


1.8 ±0.3 


4±2 


1-m 


9.0t5 


1.4 ±0.2 


76!ii 


W30H.49 


V 


443.434 


-1710.760 


15.71 ±0.04 


-42.24 


0.26 


- 


- 






- 




W3OH.50 


V 


450.607 


-1695.734 


0.62 ± 0.02 


-41.72 


0.27 


- 


- 


- 


- 


- 


- 


W30H.51 


V 


451.835 


-1673.563 


4.46 ± 0.04 


-42.24 


0.22 


1.5 ±0.2 


22 ±5 


0.55-» 


8.8- 


- 


79+" 



Notes. *"' The reference position is 0-2000 = 02''27™031833 and (52ooo = +61°52'25'.'288. **' The best-fitting results obtained by using a model 
based on the radiative transfer theory of methanol masers for F -f F^ = 1 (Vlemmings et al. 2010, Surcis et al. 201 lb). The errors were determined 
by analyzing the full probability distribution function. *^' The Zeeman- splittings are determined from the cross-correlation between the RR and 
LL spectra. *"'' The angle between the magnetic field and the maser propagation direction is determined by using the observed P\ and the fitted 
emerging brightness temperature. The errors were determined by analyzing the full probability distribution function. *''' Because of the degree of 
the saturation of these CH3OH masers TbAfl is underestimated, AVi and Q are overestimated. 



